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FOREWORD 


This document presents the results of an analysis performed 
by personnel of the Lockheed- Huntsville Research & Engineering 
Center for the Aerodynamic Systems Analysis Section of the Johnson 
Space Center, Houston, Texas, under Contract NAS9- 145 17 . The 
work was performed in support of an analysis of the Space Shuttle 
solid rocket motor exhaust plumes. 
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SUMMARY 


A parametric analysis was conducted to assess the sensitivity of the 
initial plume expansion angle of analytical solid rocket motor flow fields to 
various analytical input parameters and operating conditions. The results 
of the analysis are presented and conclusions reached regarding the sensi- 
tivity of the initial plume expansion angle to each parameter investigated. 
Operating conditions parametrically varied were chamber pressure, nozzle 
inlet angle, nozzle throat radius of curvature ratio and propellant particle 
loading. Empirical particle parameters investigated were mean size, local 
drag coefficient and local heat transfer coefficient. Sensitivity of the initial 
plume expansion angle to gas thermochemistry model and local drag coefficient 
model assumptions were determined. Initial plume expansion angle was most 
sensitive to gas thermochemistry model, propellant particle loading, and 
mean particle size assumption. Initial plume expansion angle was minimally 
sensitive to chamber pressure, nozzle throat radius of curvature ratio, 
nozzle inlet angle, local drag coefficient and local heat transfer coefficient. 

The local drag coefficient model assumption had no effect on the initial plume 
exp'ansion angle. 
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NOMENCLATURE 

Description 


nozzle area ratio 
particle drag coefficient 

particle specific heat at constant pressure 
for liquid phase at melting temperature 

particle specific heat at constant pressure 
for solid phase at melting temperature 

particle enthalpy for liquid phase at melting 
temperature 

particle enthalpy for solid phase at melting 
temperature 

Mach number 

molecular weight 

particle heat transfer coefficient 

pressure 

particle radius 

ratio of throat radius of curvature to nozzle 
throat radius 

gas static temperature 

ratio of particle mass flow rate to gas mass 
flow rate 


ratio of specific heats 
initial plume expansion angle 

nozzle lip angle 

nozzle inlet angle 


chamber condition 
exit plane condition 
ambient freestream 
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Section 1 
INTRODUCTION 


The exhaust plumes of the Space Shuttle solid rocket motors (SRMs) 
can have a significant effect (Ref. 1) on the base pressure and base drag of 
the Shuttle vehicle. Previous studies (Ref. 2) have shown that base pressure 
can be correlated to the initial plume expansion angle, 6j. These same studies 
indicate that small changes in 6. can have a significant effect on the vehicle 
base pressure. Therefore, it is necessary to predict the value of 6. as 
accurately as possible before a realistic assessment of the exhaust plume 
effect on the vehicle base pressure can be ascertained. 

Prediction of 6j requires that the nozzle exit conditions be known. This 
requires a definition of the nozzle flow field which is a function c£ the nozzle 
geometry, chamber operating conditions and particulate behavior. Nozzle 
geometry is described by analytic functions. The chamber operating con- 
ditions are obtained from analytical combustion models which utilize the 
propellant formulation and combustion pressure. However; the particle 
data (size, drag, heat transfer) are described empirically. Since the plume 
initial expansion angle has been shown to be an important correlation parameter 
for vehicle base drag, relevant questions are how does the empiricism affect 
the flow field and how sensivitive is 6j to the input data? These questions were 
the subject of an investigation which parametrically examined the sensitivity 
various input data. Input parameters examined included: 

1. Chamber Pressure 

2. Flowfield Chemistry Assumption 

3. Particle Loading 

4. Particle Size 

5. Particle Drag Coefficient 
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6. Particle Heat Transfer Coefficient, and 

7. Nozzle Geometric Modifications. 

The discussion is begun with a description of the nominal set of input 
and operating parameters for the nozzle solution. 

This report describes the results of this investigation. 
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Section 2 
DISCUSSION 

2.1 NOZZLE FLOWFIELD ANALYSIS 

The Space Shuttle SRM utilizes a solid propellant which contains 16% 
by weight of aluminum. The presence of the aluminum results in the forma- 
tion of aluminum oxide (Ai^O^) p>articulates in the exhaust. For the heavily 
aluminum loaded SRM propellant the two-phase effects on the gaseous ex- 
pansion process are significant. Therefore, it is important that the analytic 
solution of the nozzle and plume flow fields consider two- phase effects. The 
analysis of the current SRM nozzle flow field was performed using the Lockheed- 
Huntsville RAMP Two- Phase Flow program (Ref. 3). The RAMP code allows 
momentum and energy exchange between the particulate and gaseous phases, 
thus allowing the particles to affect the nozzle and exhaust plume expansion. 

The supersonic RAMP solution was initiated using a start-line generated by 
Kliegel's transonic program (Ref. 4). 

The sensitivity study was conducted by parametrically perturbing 
various input parameters to the flowfield solution about a nominal set of 
conditions. Nominal conditions were chosen to be: 

• Nozzle area ratio, A/A =7.16 

• Nozzle throat radius of curvature, R^/R^ = 2.0 

• Nozzle throat inlet angle, 0^ = 30 deg 

• Nozzle lip half angle, 0^^^ = 11.202 deg 

• Nozzle throat radius, R^ = 2.2679 feet 

• Chamber pressure, P^ = 700 psia 

• Constant gaseous thermodynamic properties 

y = 1.25 

M = 20.245 
w 
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• Particles 

loading, d) /d) = 0.4 
P 8 

mean particle size, r^ = 6.0 microns 

drag coefficient - Kliegel (Ref. 8) 

heat transfer coefficient = DraJce (Ref. 10) 


Cp^ = 0.3395 Btu/lb m-°R 
Cp = 0.2676 Btu/lb m-°R 

H = 1612.1 Btu/lb m 

Hp = 1112.8 Btu/lb m 


The nozzle flow field for these sets of nominal values provided a reference 
set of exit properties and consequently 5. as a function of expansion pressure 
ratio. The sensitivity of 6. to change in the input data was then obtained by 
parametrically varying the input data about the nominal conditions and exam- 
ming the change in 6j. 


Frozen and equilibrium gas thermochemistry models were used as 
comparative models. Data for these models were obtained using a version 
of TRAN72 (Ref. 5) which has been modified (Ref. 6) to meet the requirements 
of the RAMP Two-Phase Flow program. The taoles of thermodynamic and 
transport properties were constructed such that variations in gas properties 
due to changes in total enthalpy, entropy and temperature are considered. 


An important paramever for a two-phase calculation is the particle 
size distribution. For a nominal SRM condition a mass mean radius of 6 
microns was used. This was obtained based on a mean diameter versus 
nozzle throat diameter correlation by Delaney (Ref. 7). For comparison, 
calculations were also performed using the particle size distribution in 
Table 1. The ratio of particle mass to total propellant mass was nominally 
selected as 0.4 


For the purposes of this document, a variation in 6^ of less than one 
de^^ree is considered minimal and a variation in 6- greater than two degrees 

J 


4 


LOCKHEED • HUNTSVILLE RESEARCH A ENGINEERING CENTER 


LMSC-HREC TM D496636 


is considered signific3.nt. With the current Space Shuttle configuration, 
a one degree change in 6^ results in a change of 100 lbs in allowable payload 
weight (Ref. 7). 

2.2 CHAMBER PRESSURE EFFECT 

To investigate the effect of chamber pressure on the initial plume 
expansion angle, 6^, RAMP nozzle calculations were (generated for chamber 
pressures of 500, 600, 700, 800 and 900 psia. A figure was constructed 
with chamber to ambient pressure ratio plotted as a function of 6. for each 
chamber pressure. In Fig. 1 it is evident that chamber pressure has a ndni- 
mal effect on 6^ with a maximum variation of approximately one degree. 

2.3 EFFECT OF GAS THERMOCHEMISTRY MODEL 

The RAMP code has the capability to use several gas thermodynamic 
and transport property models. The nominal model chosen for this analysis 
is one in which the thermodynamic (y, molecular weight) and tiansport proper- 
ties (except viscosity) are held constant. The equilibrium thermochemistry model 
uses tabulated data from the TRAN72 code and allows thermodynamic properties 
to vary with changes in total enthalpy, entropy and temperature. The frozen 
thermochemistry model uses tabulated thermodynamic properties but assumes 
that no chemical reactions occur after the flow has expanded beyond a static 
to total pressure ratio specified by the user. The frozen chemistry model 
used in this analysis assumed that chemical reactions ceased at the nozzle 
throat. Nozzle calculations were performed with each chemistry model for 
a chamber pressure of 700 psia exhausting to ambient pressures ranging 
from 15.21 to 0.07 psia. 

I” F'ig* 2, chamber to ambient pressure ratio is plotted as a function of 

each chemistry model. For chamber to ambient pressure ratios less 

than 400, the use of different chemistry models did not produce any significant 

deviations in 6.. At higher chamber to ambient pressure ratios. P /p . the 
J C 00 

chemically frozen and equilibrium chemistry models produced significantly 
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different >'alucs for 6.. For a pressure ratio, P /P , of 10,000, the 

J C 00 

difference in 6. using the different chemistry models was 4 degrees. At 

J 

high pressure ratios the constant thermodynamic property model typically 

produced a 6. approximately one degree greater than the 6. produced by the 
. J J 

chemically frozen model. Above a pressure ratio of 1000, 6. values calculated 

J 

with the constant property model deviate significantly from the 6. values 

J 

calculated with the equilibrium model. The difference in initial plume ex- 
pansion for these models at P /P = 10,000 is 3 degrees. The difference in 

C 00 

y among the three models investigated appears to be the parameter which 
produces the variations in illustrated in Fig. 2. 

Nozzle calculations were performed to investigate the effect on 6^ of 
the value of Y assumed in the constant thermodynamic property model. 
Calculations were generated using values for X of 1.250 and 1.179. It is 
evident in Fig. 3 that the effect on is significant. The difference in 6. 
for the two values of X is 1 degree at P /p = 100 and 7 degrees at P /P = 

C 00 C 00 

10,000. A lower value for X produces a higher initial expansion angle. This 
trend is confirmed in Fig. 2 where the equilibrium chemistry model had the 
lowest exit plane X of the three models and produced the highest initial ex- 
pansion angles. 

2.4 NOZZLE GEOMETRY EFFECT 

Initial plume expansion angle sensitivity to two parameters defining 
nozzle geometry was investigated. The parameters were nozzle throat radius 
of curvature ratio and nozzle inlet angle. The nozzle inlet angle directly 
effects the structure of the subsonic and transonic portions of the nozzle flow 
field and the supersonic startline generated by the KUegel transonic code. 

The startline characteristics subsequently impact the supersonic expansion. 

Nozzle radius of curvature ratio and inlet angle were varied since the 
Kliegel transonic solution will not handle radius of curvature ratios below 
about 1.5 nor inlet angles greater than 45 deg. Therefore, to run two-phase 
solutions for nozzles whose geometries do not fall within the Kliegel capa- 
bility the throat geometry must be modified. For this reason nozzle geo- 
metric effects on 6. were examined. 
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Nozzle configurations with nozzle throat radius of curvature ratios 
of 1.5, 1.75, 2.0, 2.25 and 2.50 were input to the RAMP code and nozzle 
flow fields calculated with a chamber pressure of 700 psia exhausting to 
various ambient pressures. A figure was constructed with chamber to 
ambient pressure ratio plotted as a function of 6. for each throat radius of 

kl 

curvature ratio. In Fig. 4 the variation of 6. with throat radius of curvature 

J 

ratio is 0.5 degrees or less over the entire pressure ratio range. Throat 
radius of curvature ratio has a minimal effect on the initial plume expansion. 

Figures 5 and 6 are figures depicting the deviation 5^ with throat 
inlet angle for a range of operating pressure ratios.. The data pre- 
sented in Fig. 5 were calculated using one mean particle size (6 micron 
radius). The data in Fig. 6 were calculated using the particle size distribution 
of Table 1. Variation of the throat inlet angle from 20 to 45 degrees resulted 
in a 6 • variation of 0.5 to 1.0 degree at the higher pressure ratios and less 

j 

than 0.5 degree at pressure ratios less than 100. These data are presented 
in the nomographs of Figs. 5 and 6. For the calculations using one mean 
particle size, 6. was highest at a given pressure ratio for an inlet angle of 
20 degrees. For the calculations using a particle size distribution, 6^ was 
highest at a given pressure ratio for an inlet angle of 30 degrees. The vari- 
ation of 6. with throat inlet angle was considered minimal. 

2.5 PROPELLANT LOADING EFFECT 

The weight percent of solid particles in a propellant has significant 
effects on the nozzle and exhaust plume expansion of a solid propellant 
motor. To assess the sensitivity of 6^ to particle propellant loading 
CO ). nozzle flow fields were calculated for a range of propellant loadings 
(d)p/cb = 0.3 to 0.5) and chamber to ambient pressure ratios. The parametric 
calculations were generated using two chemistry models, the constant thermo- 
dynamic property model and the chemically frozen model. The figure in 
Fig. 7 indicates that at low pressure ratios the variation of 6j with 
less than 0.5 degree. However, as the pressure ratio is increased the 
variation of 6. with cb /o) increases to a significant level, e.g., a 2 degree 
variance at a pressure ratio of 6000 between calculations with a Wp/COg = 0.3 
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and calculations with a = 0.5. As cOp/Wg increases for a given pressure 

ratio, the greater the initial plume expansion. 

Results are presented in Fig. 8 for calculations using the chemically 
frozen thermochemistry model. Variations in 6j with ir^creased with 

increasing chamber to ambient pressure ratio. For a given pressure ratio 
the initial plume expansion angle increases with increasing cip/Wg. This 
effect is the result of the coupling between the gaseous and particulate phases 
which results in a higher exit pressure as particle propellant loading is 
increased. The higher exit pressure requires a greater initial plume ex- 
pansion to expand to the same ambient pressure. 


2.6 EFFECT OF MEAN PARTICLE SIZE 

An important input parameter to the analytical RAMP flowfield calcu- 
lations is mean particle size. The particle drag and heat transfer coefficients 
vary with the square of the spherical particle radius. Thus, the temperature 
and velocity lags between the gas and particle phases is a strong function of 
particle size. To investigate the deviation of 6. with different values for this 
empirically determined parameter, nozzle analyses were generated with 
various mean particle sizes. Comparative calculations were also made to 
investigate the effect of assuming one mean particle size as opposed to 
assuming a distribution of particle sizes about the given mean size. 


In Fig. 9, chamber to ambient pressure ratio is plotted as a function 
of 6. for a calculation assuming a mean particle size versus a calculation 
assuming a distribution of particle sizes (Table 1). There was no difference 
in 6. for a given pressure ratio between calculations using the different 
particle size assumptions. 

The results of nozzle calculations using different mean particle sizes 
(r =4-8 microns) are presented in Fig. 10. At a given pressure ratio, 
th*e initial plume expansion angle decreased for increasing mean particle 
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size. The difference 6. for r =4 and for r =8 was approximately 2 
degrees across the entire range of pressure ratios. Thus, there is a signi- 
ficant variation in 6j with the empirically determined mean particle size. 


2.7 EFFECT OF MEAN PARTICLE SIZE FOR DIFFERENT VALUES 
OF LOCAL DRAG COEFFICIENT 


Two empirically derived constants which are impacted by mean particle 
size and have a significant effect on the coupling between the gas and particle 
phases are the heat transfer coefficient and the local drag coefficient. An 
analysis was performed to determine the sensitivity of to different mean 
particle sizes and different heat transfer and drag coefficients. The nominal 
local drag coefficient used in calculating the results (presented in Fig. 10) 
was doubled and halved and used in nozzle calculations in which mean particle 
size and pressure ratio were parametrically varied. The results are pre- 
sented in Figs. 11 and 12, respectively. Figure 11 represents the results 
obtained with a drag coefficient double the nominal local drag coefficient. 

For chamber to ambient pressure ratios less than 700, results were similar 
to those obtained with the nominal value of drag coefficient. The intial ex- 
pansion angle decreased with increasing mean particle size at a given pressure 
ratio. However, above a pressure ratio of 700 6 . increased with increasing 

mean particle size. Differences in 6 . values ranged from 0 to 2 degrees at 

J 

pressure ratios below 700 and from 0 to 1 degree above 700. Figure 12 
represents the results obtained with a drag coefficient of one half the nomi- 
nal drag coefficient. For all values of chamber to ambient pressure ratio, 
results were qualitatively identical to those obtained with the nominal value 
of drag coefficient. The initial expansion angle decreased with increasing 
mean particle size at a given pressure ratio. Quantitatively, the variation 
of 6. with mean particle size was greater for calculations using a value of 
one half the nominal drag coefficient than for calculations using the nominal 
value. For a deviation in mean particle radius from r = 4 to r =8, 
the maximum variation in 6. was approximately 3 degrees. 
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A value of local drag coefficient of twice the nominal value resulted in 
quantitatively smaller variations in 6j than calculations using tlie nominal 
drag coefficient. Calculations using a drag coefficient one half the nominal 
value resulted in quantitatively larger variations in 6j with mean particle 
size than calculations using the nominal drag coefficient. 


2.8 EFFECT OF MEAN PARTICLE SIZE FOR DIFFERENT VALUES 
OF LOCAL HEAT TRANSFER COEFFICIENT 


The nominal local heat transfer coefficient used in calculating the 
results presented in Fig. 10 was doubled and halved and used in nozzle 
calculations in which mean particle size and pressure ratio were parametri- 
cally varied. The results are contained in the nomographs of Figs. 13 and 
14, respectively. Figure 13 represents the results obtained with a local 
heat transfer coefficient double the value of the nominal heat transfer 
coefficient. For chamber to ambient pressure ratios greater than 700, 
results were similar to those obtained with the nominal value of heat trans- 
fer coefficient. The initial expansion angle decreased with increasing mean 
particle size at a given pressure ratio. The variation in 6j with mean particle 
size was negligible at a chamber to ambient pressure ratio of 700. However, 

below a pressure ratio of 700 6. increased with increasing mean particle 

J 

size. Differences in 6. values for various mean particle sizes varied from 
0 to 0.5 degrees across the entire pressure ratio range. 


Figure 14 represents the results obtained with a local heat transfer 
coefficient of one-half the value of the nominal heat transfer coefficient. 

For all values of chamber to ambient pressure ratio, results were qualita- 
tively identical to those obtained with the nominal value o^ heat transfer 
coefficient. The initial expansion angle decreased with increasing mean 

particle size at a given pressure ratio. Quantitatively, the variation in 6. 

J 

with mean particle size was greater for calculations using the value of one- 

half the nominal heat transfer coefficient than the calculations using the 

nominal heat transfer coefficient. For a variation in mean particle radius 

from r = 4 to r =8, the maximum deviation in 6. was approximately 2 
mp mp j ^ 
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do^rcH‘S for calculations using values of one-half the nominal heat trans- 
fer coefficient. 

2.9 EFFECT OF DRAG COEFFICIENT MODEL 

Several empirical models are currently in use for calculating particle 
drag coefficients. The models were developed by Kliegel (Ref. 8) and Crowe 
(Ref. 9), respectively. An analysis was conducted to assess the sensitivity 
of 6. to different drag coefficient models. In Fig. 15 chamber to ambient 
pressure ratio is plotted as a function of 6. for calculations using drag co- 
efficient models developed by Kliegel and Crowe. There was no difference 
in 6. at a given pressure ratio for calculations using different drag coefficient 
models. These calculations were performed assuming a mean particle size 

with r =6 microns. 

mp 

An analysis was conducted to investigate the sensitivity of 6^ to different 
values of drag coefficient. Using the Kliegel drag model and parametrically 
varying pressure ratio, nozzle calculations were generated using different 
values of drag coefficient (0.5 x Kliegel Cj^ to 2.0 x Kliegel Cj^). The results 
of these calculations are represented by the nomograph of Fig. 16. For a 
given pressure ratio, the initial plume expansion angle increases slightly for 
increasing values of drag coefficient. Over the range of pressure ratios 
investigated, the maximum variation in was less than one degree. The 
increased initial plume expansion angle with increased value of drag coefficient 
resulted from the coupling between the gas and particulate phases. The in- 
creased particle drag coefficient increased the momentum loss of the gas 
resulting in a lower gas velocity and higher pressure at the nozzle exit plane. 
The higher exit pressure produced a greater initial plume expansion. The 
change in 6. with the value of drag coefficient (assuming a mean particle size) 
was minimal over the range of drag coefficient values considered. 

The sensitivity of 6. to different drag coefficient models assuming a 
distribution of particle sizes was investigated. The distribution of particle 
sizes in Table 1 was used in nozzle calculations with drag coefficient models 

1 1 
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by Kliegel and Crowe. In Fig. 17 chamber to ambient pressure ratio is 
plotted as a function of 6j for calculations using different drag coefficient 
models. At a given pressure ratio, there was no difference in 6j for nozzle 
calculations using different drag coefficient models. 

An analysis was conducted to investigate the sensitivity of 6j to different 
values of drag coefficient using a distribution of particle sizes. Using the 
Kliegel drag model and parametrically varying pressure ratio, nozzle cal- 
culations were generated using different values of drag coefficient (0.5 x 
Kliegel to 2.0 x Kliegel Cj^). The results of these calculations are pre- 
sented in Fig. 18. Across the range of pressure ratios investigated, the 
initial plume expansion angle increased with increasing drag coefficient. 

As pressure ratio increased, the difference in 6^ for calculations using 
various values of Cj^ increased. Quantitatively, for the range of drag 
coefficients used, the maximum deviation in 6^ at given pressure ratio 
was approximately one degree. The quantitative and qualitative results 
were the same as the results obtained from calculations which assumed 
a mean particle size. Changes in 6^ with values of local drag coefficient 
were considered minimal. 

2.10 EFFECT OF VALUE OF HEAT TRANSFER COEFFICIENT 

Nozzle calculations were generated to assess the effect of the value 
of heat transfer coefficient on the initial plume expansion angle. Pressure 
ratio was parametrically varied for calculations using values of the local 
heat transfer coefficient ranging from 0.5 x Drake Q to 2.0 x Drake In 
Fig. 19, the variation of 6. with heat transfer coefficient at a given pressure 
ratio is illustrated. The initial plume expansion angle increases with in- 
creasing local heat transfer coefficient for a given chamber to ambient 
pressure ratio. For the range of pressure ratios and heat transfer coeffi- 
cients investigated, difference in 6^ at a given pressure ratio ranged from 
0.5 to 1.0 degree. 
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To investigate the combined effect of heat transfer coefficient and 
thermochemistry model, parametric calculations were performed assum- 
ing the flow to be chemically frozen at the nozzle throat. The value of the 
local heat transfer coefficient was varied from 0.5 x Drake 0 to 2.0 x 
Drake Q. The results of the calculations are presented in Fig. 20. For 
a given pressure ratio. 6^ increases with increasing value of heat trans- 
fer coefficient. The maximum deviation in 6. at a given p ressure ratio 
was 2.5 degrees. Typically, the variation in 6. was greater in the calcu- 
lations using the chemically frozen thermochemical model than in the 
nominal calculations assuming constant thermodynamic properties. 


To investigate the combined effect of heat transfer coefficient and parti- 
cle size distribution, parametric calculations were performed assuming the 
particle size distribution of Table 1 rather than the nominal mean particle 
radius of 6 microns. As in the previous cases, the value of the local heat 
transfer coefficient was varied from 0.5 x Drake Q to 2.0 x Drake Q. 
Qualitatively and quantitatively the results presented in Fig. 21 are similar 
to those in Fig, 19 obtained with a mean particle radius of 6 microns. At 
a given chamber to ambient pressure ratio, 6^ increased with increasing 
heat transfer coefficient. The variation of 6^ with heat transfer coefficient 
over the range of pressure ratios was approximately 0.5 to 1.0 degree. It 
was concluded that the differences in particle size distributions investigated 
did not affect the variation of 6. with heat transfer coefficient. 

J 


2.11 


COMBINED EFFECT OF PARTICLE 
TRANSFER COEFFICIENT 


DRAG AND HEAT 


The sensitivity of 6j to simultaneous changes in local drag and heat 
transfer coefficients was investigated. Nozzle calculations were performed 
with; (1) a local drag coefficient of twice the Kliegel drag coefficient and a 
local heat transfer coefficient of twice the Drake heat transfer coefficient 
(case 1), (2) a nominal heat transfer coefficient and a drag coefficient of 
one-half the Kliegel coefficient (case 2), (3) a nominal drag coefficient and 
a heat transfer coefficient of one-half the Drake coefficient (case 3). 
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Results of these calculations are plotted comparatively in Figs. 22 through 
25. In Fig. 22 the parametric variation of 6. with pressure ratio for cases 
1 and 3 is plotted along with the results of calculations using nominal values. 
Over most of the pressure ratio range, the nominal curve fell between the 
curves for calculations using one-half the Drake heat transfer coefficient 
and for calculations using twice the value of the Kliegel drag and Drake heat 
transfer coefficients. The maximum difference in 6^ for the two later curves 
was approximately two degrees. The results of calculations for cases 1 and 2 
and nominal conditions are presented in Fig. 23. There was no difference in 
values of 6. at a given pressure ratio for the nominal curve and the curve 

J 

generated with one-half the value of the local drag coefficient. This result 
is verified by results presented in Fig. 16. At pressure ratios below 500, 
the curves for all three cases coincided. Above a pressure ratio of 500, 
calculations using drag and heat transfer coefficients twice the nominal values 
diverged from the nominal curve. At a given pressure ratio above 500, the 
initial plume expansion was greater for calculations using twice the nominal 
drag and heat transfer coefficients. This trend agrees qualitatively with the 
results presented in Figs. 16 and 19 in which the drag and heat transfer co- 
efficients, respectively, were increased separately. 

Figure 24 compares the results of three calculations: (1) using a nomi- 

nal value of heat transfer coefficient and a drag coefficient of twice the nomi- 
nal value; (2) using twice the nominal values of heat and drag coefficients; and 
(3) using nominal values. The curves of chamber-to-ambient pressure ratio 
as a function of 6. for all three conditions coincide at pressure ratios below 

ij 

500. Above a pressure ratio of 500, the curves for the off-nominal conditions 
diverged from the curve generated with nominal values for drag and heat 
transfer coefficients. Both curves generated with twice the nominal value 
of drag coefficient coincided over the entire pressure ratio range investigated. 
Thus the drag coefficient appeared to be a stronger driving function of 6j than 
the heat transfer coefficient. Qualitatively, the trends are similar to the 
results presented in Figs. 16 and 19 in which heat transfer and drag coefficients 
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were increased and produced larger 6j values at a given pressure ratio. 
Above a pressure ratio of 500, calculations using values of drag and heat 
^*"^^sfer coefficients of twice the nominal values produced 6. values of up 
to 1 degree greater at a given pressure ratio than calculations using the 
*^orni*^3’l values. The maximum difference in 6j (1 degree) between nominal 
and off-nominal calculations (assuming values of heat transfer and drag 
coefficients twice the nominal values) was the same order of magnitude as 
the difference in 6. between nominal calculations and calculations generated 

J 

with the coefficients doubled separately (Figs. 16 and 19). 
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Section 3 
CONCLUSIONS 

The sensitivity of the initial plume expansion to various operating and 
input parameters was investigated. The variation of 6. with changes in 
chamber pressure at constant pressure ratio was minimal. At low chamber- 
to ambient pressure ratios, the use of different chemistry models did not 
produce any variation in 6^. The difference in X among the three chemistry 
models resulted in large difference in 6^ at the high pressure ratios. The 
largest difference in 6^ at a given pressure ratio occurred between calcu- 
lations using a constant property model and calculations using an equilibrium 
chemistry model. Large differences in 6^ (7 degrees) were obtained when 
different values of X were used in calculations with a constant property 
chemistry model. At a given pressure ratio, as X decreases the initial 
expansion angle, 6^, increases. Changes in nozzle throat radius of curvature 
ratio and nozzle inlet produced minimal changes in 6.. At pressure ratios 
above 1000, propellant particle loading produced significant differences in 
6j. The initial plume expansion angle increased with increasing propellant 
loading at a given pressure ratio. There was no difference in 6. between 
calculations assuming a mean particle size and calculations assuming a 
particle size distribution. As mean particle size was varied from r =4 
’’mp “ initial plume expansion angle decreased by Z degrees. ^The 

assumed mean particle size thus has a significant effect on 6.. There was 
no difference in 6^ at a given pressure ratio for calculations using different 
drag coefficient models, e.g., Kliegel and Crowe. Parametric variation of 
local drag coefficient from 0.5 x Kliegel to 2.0 x Kliegel Cj^ produced a 
maximum 6^ variation of less than one degree. The parametric variation of 
heat transfer coefficient indicated increasing 6j with increasing values of 
local heat transfer coefficient. The deviation of 6. with heat transfer co- 
efficient was considered minimal. The effect on initial plume exfjansion 
angle from combined changes in local drag and heat transfer coefficients 
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was investigated. Qualitatively and quantitatively, the results were similar 
to those obtained in calculations in which only one of the coefficients was 
changed at a time. It was concluded that the coupling between the drag and 
heat transfer calculations was minimal. 

In summary the following qualitative effects on initial plume expansion 
were observed for each parameter: 

• Chamber pressure — minimal effect 

• Chemistry model — large variation in 6^ 

• Nozzle throat radius of curvature ratio and throat 
inlet angle — minimal effect 

• Mean particle size assumption versus particle dis- 
tribution assumption — no effect 

• Mean {larticle size — large variations in 

• Propellant particle loading (ob 
variations in 6. 

J 

• Drag coefficient model (Kliegel versus Crowe) — 
no effect 

• Drag coefficient value — minimal effect 

• Heat transfer coefficient value — minimal effect 

• Combined changes in value of drag and heat 
transfer coefficients — minimal effect. 


p/(bg) — significant 
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Variation of 6j with Change in Nozzle Throat Radius of Curvature Ratio 
and Ambient Pressure Ratio 
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Fig. 19 - Variation of 6. with Change in Local Heat Transfe 
Coefficient and Ambient Pressure Ratio 


LOCKHEED- HUNTSVILLE RESEARCH & ENGINEERING CENTER 

































3303.5 


Fig. 20 - Variation of 6j with Change in Heat Transfer Coefficient and 

Ambient Pressure Ratio Assuming Chemically F rozen Flow 
from the Nozzle Throat 
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